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ABSTRACT 



We present the analysis of simultaneous multicolour 
uvylc photometry and low-resolution spectroscopy for 
the rapidly rotating (3 Cephei star S Y Equ. From the pho- 
tometric time series, we confirm the dominant pulsation 
frequency, /i = 6.029 d _1 , and we find an evidence for 
two additional modes. In spectroscopy, the highest peak 
occurs at / a = 0.197 d _1 or its alias 0.803 d -1 . It can 
be interpreted either in terms of a binary motion or as 
the g-mode pulsation. In addition, we reveal the pulsa- 
tion mode with frequency of about 6.029 d _1 , i.e. the 
same which dominates photometric variations, and a few 
new candidates. For the dominant frequency we obtain 
mode identification from the combined photometric and 
spectroscopic observations. From non-adiabatic pulsa- 
tion calculations, we show that the frequency of the domi- 
nant mode in S Y Equ is consistent with the stellar models 
of much lower effective temperatures than used in many 
papers. 
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(3 Cephei-type variables; pulsation; rotation. 



1. INTRODUCTION 



SY Equ (HD 203664) is the early-B spectral type star, 
B0. 5 Illn, classified as the (3 Cephei variable by from 
the Hipparcos photometry. The star is one of the fastest 
rotators among (3 Cephei pulsators, having projected ro- 
tational velocity of about 200 km s _1 0. Moreover, it 
is a high Galactic latitude star with b = —21 .5° which lo- 
cates it more than 1 kpc out of the Galactic plane. 

There are several determinations of stellar parameters for 
SY Equ. derived log T cS = 4.447 dex, log g = 3.7 dex 
from the high-resolution optical spectra and estimated a 
mass of 14 Mq from the evolutionary tracks of On 
the other hand, using Geneva photometry, determined 
logTeff = 4.47 ± 0.01 dex and logg = 3.9 ± 0.3 dex. 

* Based on photometry obtained at the Bialkow Observatory (Uni- 
versity of Wroclaw) and on spectroscopy obtained at the David Dunlap 
Observatory (University of Toronto) 



Much lower value of effective temperature, logTeff = 
4.401 dex, was obtained by |6] from Stromgren photom- 
etry and equivalent width of the He line. We adopt here 
stellar parameters derived from the International Ultravi- 
olet Explorer (IUE) spectra and visual photometry by 0: 
logT eff = 4.388 ± 0.035 and \ogg = 3.905 ±0.1. These 
values result in a mass estimate of 1 1 .5 ± 1.5 Mq if the 
metallicity Z = 0.02 is assumed. The metallicity in terms 
of [m/H], obtained from the IUE/INES data, amounts to 
0.00 ± 0.21 dex, which corresponds to Z = 0.02^°;°^. 
We prefer to rely on this determination because it uses ul- 
traviolet spectral region, where B-type stars emit most of 
their energy. 

In this note we give only a brief outline of our results. 
The full analysis will be published elsewhere. 



2. OBSERVATIONS 



2.1. Photometry 

The new CCD photometry of SY Equ was carried out in 
Bialkow Observatory (Wroclaw University) during five 
nights between August 14 and September 3, 2004. We 
used Stromgren uvy and Cousins Iq passbands. The ob- 
tained frames were calibrated in a standard way and then 
reduced using the profile-fitting and aperture photometry 
of the Daophot package Jl2ll . 



2.2. Spectroscopy 

Spectroscopic observations of SY Equ were obtained 
with the Cassegrain spectrograph attached to the 1.88-m 
telescope of the David Dunlap Observatory (University 
of Toronto). The data consist of 171 spectra taken on 12 
nights between August 8 and September 1, 2004, i.e., si- 
multaneously with the Bialkow photomety. The spectra 
were sampled over the range 4000^4600 A with the 600 
lines per mm grating giving a dispersion of 0.62 A per 
pixel. The exposure times were 15 minutes long. The re- 
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Figure 1. The mean spectrum of SY Equ in the observed 
range obtained by averaging of 171 normalized individ- 
ual spectra. IB stands for interstellar band. 



ductions were made with standard IRAF 1 routines, which 
include cosmic ray removal, bias correction, flat-fielding, 
and wavelength calibration. The spectra were moved to 
the heliocentric frame. In Fig.^ we show the mean spec- 
trum of SY Equ. 



3. FREQUENCY ANALYSIS 



The photometric light curve of SY Equ is dominated 
by the mode with frequency / = 6.029 d _1 . The 2004 
Bialkow vbylc photometry folded with the period of 
0.165871d= 1/fis shown in Fig. |2] In order to search for 
possible low-amplitude modes, we combined three pho- 
tometric datasets: Geneva F-band data of 01, ASAS-3 
y-band photometry Jill EU and Bialkow 2004 data in 
the Stromgren y band. These data cover the time interval 
of almost 9 years, between 1997 and 2005. In Fig. [3] we 
plotted Fourier periodograms of the combined Vly data 
showing consecutive steps of prewhitening. We see that 
except for the well-known dominating mode, at least two 
low-amplitude modes are present. Their frequencies are 
equal to 8.360 and 7.821 d _1 , but due to the severe alias- 
ing, their values are not certain. 

For the spectroscopic data, we searched for intrinsic vari- 
ations in the first four moments of four spectral lines: 
H 7 , H<5, He 4471, and He 4387, labeled in Fig.[T] More- 
over, we did the same for radial velocities obtained by 
means of the cross-correlation function (CCF) technique 
using the whole observed spectrum excluding the region 
of the intestellar band (IB, Fig. 0. The frequencies of 
the terms we found are summarized in Tab. 1. In all 
lines and the CCF radial velocity, we have found low- 
frequency variation at 0.197 d — 1 > which did not appear 
in photometry, and the peak known already from pho- 
tometry, ^6.029 d _1 . The radial velocity amplitudes of 

1 IRAF is distributed by the National Optical Astronomy Observato- 
ries, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Sci- 
ence Foundation. 
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Figure 2. Phase diagrams of the Bialkow 2004 SY Equ 
light-curves folded with the period of 0.165871 d. Going 
from top to bottom, data in u, v, y, and Iq photometric 
passbands are shown. 
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Figure 3. Fourier periodogram of the of the combined V/y 
photometric data of SY Equ (top) and after consecutive 
steps of prewhitening (lower panels). 



these two frequencies are 15 and 12 km s _1 , respec- 
tively. Additionally, we found low-amplitude changes 
with 4.42 d _1 and 0.61 d -1 . The latter could be the 
rotation frequency of SY Equ. In Table 1 we give fre- 
quencies which were identified after subsequent stages 
of prewhitening for three moments of two lines (H7 and 
H<5) and y rac i(CCF). The important results is that we did 
not find the variation with frequency 2/i, where f\ is the 
frequency of the dominant mode, in the second moments 
of the analyzed four lines. This leads immediately to the 



Table 1. Frequencies [in d~ ] detected with signal-to-noise (S/N) ratio greater than 4 in the first three moments of 
and H8 lines, and CCF radial velocities. The frequencies are given in the order of detection. Numbers in parentheses 
denote the r.m.s. errors of the last digits. 
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V Tad (CCF) 


6.0318(39) 


0.1970(15) 
6.0267(16) 
0.6094(21) 
9.0304(28) 


0.1921(27) 
6.0283(30) 


0.4683(42) 
6.0284(37) 


0.1970(16) 
6.0253(18) 
4.4205(33) 


6.0354(14) 
0.3532(21) 


6.0241(12) 
0.1978(19) 
4.4281(28) 
0.6026(36) 



conclusion that the dominant frequency has to be a non- 
axisymmetric mode (m ^0). 

Recently 0] suggested two additional modes with fre- 
quencies of 6.82902 and 4.81543 d _1 , from the analysis 
of the Geneva photometry. The former is probably an 
alias of our mode with 7.821 d _1 , the presence of the lat- 
ter was not confirmed in our analysis. 



4. MODE IDENTIFICATION 



We determined the amplitudes and phases of the light 
curves in uvylc bands and of the radial velocity varia- 
tions. Using these data for the dominant frequency, we 
identified its spherical harmonic degree, t, by means of 
the method invented by Daszyhska, Dziembowski and 
Pamyatnykh 0H|. The method combines photometry 
and spectroscopy and consists of simultaneous determi- 
nation of the £ value and the pulsation nonadiabatic pa- 
rameter /. The / parameter describes the ratio of the 
relative luminosity variation to the radial displacement of 
the surface, and it gives information on subphotospheric 
layer. In the case of 8 Scuti stars, this new diagnostic 
tool probes the efficiency of convection |4], while in the 
case of (3 Cep stars, it provides a stringent probe of stellar 
opacities 01 . From this method, we get the unambiguous 
£ = 2 identification. In Fig. |4] we show the \ 2 value 
as a function of I. The horizontal line at x 2 =2.1 corre- 
sponds to 95% confidence level. In the next step, we iden- 
tified the azimuthal order, m, from spectroscopy, using 
monochromatic amplitude and phase diagrams across the 
line profiles. In Fig.|5]we plot these diagrams across H7. 
As one can see, the amplitude reaches a very sharp max- 
imum in the line center, which indicates \m\ = 2 as the 
most probable identification. The sign of m is obtained 
from the phase changes. Because the phase decreases 
across the line profile with increasing wavelengths, the 
mode must be prograde, m — +2. Thus, the dominant 
frequency of SY Equ is a sectoral prograde mode with 
1 = 2. Such mode identification excludes low inclination 
angles because sectoral modes are poorly seen when the 
star is observed near pole-on. 




Figure 4. The value of x as a function of £ obtained 
for five models of SY Equ. The models are located in the 
error box defined in Sect. 1. 



5. PULSATION MODELING 



Having identified the geometry of the main pulsation 
mode in terms of £ and m parameters, we may compare 
the observed frequency of this mode with the theoretical 
predictions. For this purpose, stellar evolutionary models 
were calculated using Warsaw-New Jersey evolutionary 
code, assuming the OPAL opacities and the standard solar 
mixture of elements. Linear nonadiabatic pulsation cal- 
culations were performed with the use of Dziembowski's 
code, which includes the rotation effects on oscillation 
frequencies up to the second order in the rotational ve- 
locity. Although SY Equ rotates very fast (at least 35% 
of the break-up velocity), we can neglect deviations from 
the spherical symmetry due to effects of centrifugal force, 
because the ratio of the rotation frequency to pulsation 
frequency is still much less than 0.5. 

In Fig.|6j unstable modes with £ = 0, 1, and 2 for models 
with a mass of 12 M , evolving from ZAMS to TAMS, 
are shown. The initial hydrogen abundance was X = 0.7 
and the heavy element abundance Z = 0.02. We con- 
sider two values of the rotational velocity, v ro t= and 
200 km s _1 . The vertical lines mark the range of effective 
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Figure 5. The monochromatic amplitude ( top) and phase 
(bottom) diagrams for H^j. The line profile is plotted with 
a thin line and the amplitudes and phases with a thick 
one. Errors are shown as grey shadows. Right ordinate 
is given in flux units and refers to the line profile. 



temperatures determined from the IUE spectra. The rota- 
tional velocity of SY Equ is equal to at least 200 km s _1 , 
which gives rotational frequency of ~0.6 d _1 assuming 
the stellar radius of about 6.5 Rq. If the frequency of 
6.029 d" 1 (solid horizontal line) is a sectoral prograde 
mode, we have to compare theoretical values with the 
centroid of the £ = 2 mode which should be placed at 
about 4.8 d _1 (dotted horizontal line). As one can see 
from Fig. [6] for Z = 0.02 the pulsation instability begins 
at logT e ff of about 4.41 dex, which is much lower than 
the effective temperature of SY Equ used by some au- 
thors, e.g. 01 and |7|. Moreover, we checked the metal- 
licity and mass effects on the instability domain consider- 
ing several values of Z and M. Changing Z from 0.015 
to 0.025 shifts the hottest unstable mode from log T e g s» 
4.40 to 4.42 dex. Changing mass from 1 1 to 13 M® gives 
exactly the same effect on the hot edge of instability. 



6. CONCLUSIONS 



We presented the results of the frequency analysis of both 
the photometric and spectroscopic time series of S Y Equ. 
In both datasets, we found the dominant mode at f-y = 
6.029 d _1 , and several new candidate frequencies. More- 
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Figure 6. Evolution of unstable modes with i = 0, 1, and 
2 in a sequence of (3 Cephei star models with mass of 12 
Mq, assuming Z = 0.02. The models cover the phase 
from ZAMS to TAMS. The two panels show the effects of 
rotational velocity (0 and 200 kms~ x ) on the pulsation 
frequencies. 



over, there is a low-frequency variation at f a = 0. 197 d _1 
detected in radial velocities, which can be an evidence 
of the binary motion or the excitation of a gravity mode. 
For the dominant mode we confirmed the £ = 2 identifica- 
tion of [2] and identified the azimuthal order, m = +2. It 
means that the dominant frequency is a sectoral prograde 
mode with £ =2. However, one has to be aware that at 
V m t ~ 200 km s ~ 1 , the rotational mode coupling can take 
place, if the frequency difference between modes j and k 
is of the order of angular velocity of rotation, and if the 
spherical harmonic indices satisfy the relation lj = i^i 
2, see d. 

Recently, 0] made an attempt of comparisons of theoret- 
ical frequencies with the observed value of the dominant 
mode, but the authors totally ignored the effects of rota- 
tion and the instability condition. Our analysis showed 
that only for lower effective temperatures we enter the 
instability domain. These temperature values agree very 
well with the range determined from the ultraviolet IUE 
spectra. Moreover the rotation effects on frequencies can- 
not be neglected for such rapid rotator as S Y Equ. 
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